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We discuss the relation between entanglement and criticality in translationally invariant harmonic 
lattice systems with non-randon, finite-range interactions. We show that the criticality of the 
system as well as validity or break-down of the entanglement area law are solely determined by the 
analytic properties of the spectral function of the oscillator system, which can easily be computed. 
In particular for finite-range couplings we find a one-to-one correspondence between an area-law 
scaling of the bi-partite entanglement and a finite correlation length. This relation is strict in the 
one-dimensional case and there is strog evidence for the multi-dimensional case. We also discuss 
generalizations to couplings with infinite range. Finally, to illustrate our results, a specific ID 
example with nearest and next-nearest neighbor coupling is analyzed. 
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Due to the development of powerful tools to quantify 
entanglement there is a growing interest in the relation 
between entanglement and criticality in quantum many- 
body systems. For a variety of spin models it was shown 
that in the absence of criticality, there is a strict relation 
between the von-Neumann entropy of a compact sub-set 
of spins in the ground state and the surface area of the 
ensemble. E.g. it was shown in Q, 0, 0, 01 that the 
entanglement in non-critical one-dimensional spin chains 
approaches a constant value, while it grows logarithmi- 
cally in the critical case, where the correlation length 
diverges. Employing field theoretical methods it was ar- 
gued that in d dimensions the entropy grows as a poly- 
nomial of power d — 1 under non-critical conditions, thus 
establishing an area theorem. A similar relation was sug- 
gested for harmonic lattice models in @ and @. Very 
recently, employing methods of quantum information for 
Gaussian states, Plenio et al. [7j gave a derivation of the 
area theorem for harmonic lattice models with nearest- 
neighbor couplings. All these findings suggest a general 
correspondence between entanglement and criticality for 
non-random potentials. Yet recently special cases have 
been found for spin chains with Ising-type interactions || 
and for harmonic lattices systems 9] where the correla- 
tion length diverges but the entanglement obeys an area 
law. Thus the relation between entanglement scaling and 
criticality remains an open question. It should also be 
noted that in disorderd systems, i.e. systems with ran- 
dom couplings the relation between entanglement area 
law and criticality is broken. 

In the present paper we show that for harmonic lattice 
systems with translational invariant, non-random, and 
finite-range couplings both entanglement scaling and crit- 
icality are determined by the analytic properties of the 
so-called spectral function. For finite-range interactions 
we find that the properties of the spectral function lead to 
a one-to-one correspondence between entanglement and 
criticality. To illustrate our results we discuss a specific 



one-dimensional example with nearest and next-nearest 
couplings. Despite the finite range of the coupling this 
model undergoes a transition from area-law behavior to 
unbounded logarithmic growth of entanglement. 

Let us first consider a one-dimensional system, i.e. a 
chain of N harmonic oscillators described by canonical 
variables (qi,Pi), i = 1,2,...N. The oscillators are cou- 
pled by a translational invariant quadratic Hamiltonian 



where V is a real, non-random, symmetric matrix with 
positive eigenvalues. For a translational invariant system 
V is a Toeplitz matrix, i.e. its elements depend only on 
the difference of the indexes Vij =Vh = V-h- For a finite 
system translational invariance implies furthermore peri- 
odic boundary conditions V& = V/v-fc. We assume in the 
following that the interactions are of finite range, i.e that 
Vk = for k > R, where R is a finite number indepen- 
dent on N. As we will show at the end of the paper some 
generalizations to infinite range couplings are possible. 
Being positive definite, V has a unique positive square 
root V 1 / 2 and its inverse V~ x l 2 which completely deter- 
mine the ground-state in position (or momentum) rep- 
resentation, *o(Q) ~ (dety 1 /2)i/4e Xp {-i(Q|yV2|Q\} ) 

m 

The most important characteristics of the oscillator 
system is the spectrum of V. Since V is a ciculant matrix 
its eigenvalues can be expressed in terms of complex roots 
of unity Zj = exp{i2irj / N} — exp{i6j}, (j = 1 , . . . , N): 

a, = E v * fe) fe = \ v « + \ £ v * ( c ^') fe(2) 

k=0 fc=-(il-l) 

Eq.(J2J together with the positivity of V permits the 
representation Xj = h (zj) = \h 2 (zj)\ where h(zj) is 
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a polynomial of order (R — l)/2 in Zj (assuming for 
simplicity that R is an odd number). Thus ~ 

n[=i I- 2 — Zl\- where Z\ = exp{«a;} are the zeroth 
of h(z) which are either real or complex conjugate pairs 
with \zi\ > 1 [lo|. Let Q < R be the number of 
real zeroth z r with multiplicity m r € {0, 1, . . . }. Then 



IELlh - = ntt(2 - 2C0S(^ - Or))' 

Q 

A(zj) = Aj = Xo(zj) TT(2 — 2cos(0j — a r ) 
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and 



(3) 



\(z) is the so-called spectral function. Xq{z) is called 
the regular part of A. It is a polynomial of the complex 
variable z which has zeroth outside the unit circle. As 
a consequence its inverse Aq X (z) is analytic on and in- 
side the unit circle. — 2003(6* — a r )) m '' is called 
the singular part. If A is singular, i.e. if in the thermo- 
dynamic limit V has eigenvalues arbitrarily close to zero, 
the total Hamiltonian, eq.JQ), has a vanishing energy gap 
between the ground and first excited state. 

To evaluate sums of eigenvalues in the limit N — > 00, 
one can interpret eq.Q as a Fourier series of X(6). Thus 
Vfe = ^ J o ^d0 X(0) e~ t6k . This integral representation is 
also valid for a finite number TV of oscillators up to an 
error 0(1/N) as long as k < (N + l)/2. Due to the peri- 
odic boundary conditions V ±x / 2 are also Toeplitz matri- 
ces, and their elements V^ 1 ^ 2 — V^ 1 ^ 2 can be expressed 
in terms of A ±x / 2 for k < (N + l)/2 



V ±l/2) - 
k 



= _Lp 9 ^ (e)e - 



■iS k 



(4) 



Since for the spatial correlation of an oscillator system 
holds (qiqi+i) ~ V t 11], the analytic properties of 
A -1 / 2 determine the spatial correlation length £: 
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If some derivative of X~ 1 ^ 2 (9), say the mth one, does not 
exists, partial integrations shows that the integral has a 
contribution proportional to l~ rn . In this case the corre- 
lation length £ is infinite, defining a critical system. If 
on the other hand A -1 / 2 ((9) is smooth the integral decays 
faster than any polynomial in In this case the cor- 
relation length is finite, corresponding to a non-critical 
system. From the form of X(8) given in eq.© it is clear- 
that a regular spectral function implies a finite correla- 
tion length, i.e. a noncritical behavior and a singular one 
an infinite correlation length, i.e. a critical behavior. 

In the following we will show that the analytic prop- 
erties of A also determine the entanglement scaling of 
the oscillator system. The bi-partite entanglement of a 



compact block of Ni oscillators (inner partition X) with 
the rest (outer partition O) is determined by the N\- 
dimensional sub-matrices A and D [3, la, LD, LUj 



V 
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C and F are here {N - N x ) x (N - N\) matrices. The 
entropy is given by the eigenvalues Hi > 1 of the matrix 
product A ■ D Q: 



(7) 



where f(x) — ^4p- In ^4p- — In Despite the sim- 
plicity of its form, J7J) cannot be evaluated in general. 
This is in contrast to spin systems where A ■ D is itself a 
Toeplitz matrix 0,0]. 

An upper bound to S can be found from the logarithmic 
negativity ln||p r ||, where p T is the partial transpose of 
the total ground state p and 1 1 • 1 1 denotes the trace norm. 
As shown in 0, the logarithmic negativity is bounded 
by the square root of the maximum eigenvalue of V and 

— 1/2 

a sum of absolute values of matrix elements of V,, 
between all sites i el and j € O. 
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A lower bound to the entropy can be found making use 



of 2±1 In 2±1 - In > In x . This yields 



1 Nl 
S > 2$I ln/ii 

i=l 



-ln( detA-D 



(9) 



This estimate has a simple and very intuitive meaning. 
To see this we first note that the matrix D can be ex- 
pressed in the form D = (A - B ■ C ■ B T y 1 . Thus 



S > -ilndct^-B-Cr 1 



• In det 
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det F det D 
det V 1 / 2 
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where the last equation was obtained by expressing A 
in terms of D,E, and F. The last line of (|10|l is just 
Shannon's classical mutual information I(Qi : Q2) or 
I {Pi : Pa) respectively, where Qi = (qi, q%, . . . , q^) an d 
Q2 = (ffJVi+i) • • • j Qn) are the position vectors of the two 
subsystems and P\,i the respective momentum vectors. 
I(Qi : Q2) is defined as 



I(Qi :Qi)= / d N Qp{Q l ,Q 2 ) In 



p{Qi,Q2) 

Pl(Ql)P2(Q2 



(11) 
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where p(Qi,Q%) = \^o\ 2 is the total and ^1,2(^1,2) the 
reduced probability density in position space. Straight 
forward calculation shows 

. „ , 1 , det A- det C „ , „. 
I(Q i: Q 2 ) = -ln dety _ 1/2 <5. (12) 

In order to evaluate Shannon's mutual information in 
the form given in eq.JUJ) we want to make use of the 
asymptotic properties of Toeplitz matrices. For this we 
note that since V^ 1 / 2 are Toeplitz matrices, so are A and 
D. Their elements A k , and D k can be obtained from 
A±V2 by if N x < (N + l)/2. 

If M9) is regular, we can apply the strong Szego theo- 
rem Il2], which states: 



00 

det(D) -> exp{c iVi+^fc|c fc | 2 }, 



fc=0 



(13) 



for iVi — > 00. Here the Ck are Fourier-coefficients of 
In A" 1 / 2 (gi), i.e. c k = ± J 27I d9 In A 1 / 2 ((9) e" iefc . Noting 
that the corresponding coefficients for A have opposite 
sign, we find the lower bound 



S > - ln(det(A)) 

£1 



■ln(det(£»)) 
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To find an upper bound to S we make use of eq.lJSJ. 
For a finite-range interaction there is always a maximum 



eigenvalue Am ax . Furthermore since A (9) is smooth, 
eq.Q implies an exponential bound to the matrix ele- 



ments of V 1 / 2 . I.e. 



\Vr 1/2 \ < Xexp{-a|i-j|}, for 



• j\ < (N + l)/2, where K, a > 0. With this we find 



(N+l)/2 
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for N,Ni — > 00. Thus 5 has also a finite upper bound 
in ID. One recognizes that for one dimensional harmonic 
chains with a regular spectral function X(9) the entropy 
has a lower and an upper bound independent on the num- 
ber of oscillators, which implies an area theorem. Fur- 
thermore, as shown above, the spatial correlation length 
is finite, i.e. the system is non-critical. 

Let us now consider a singular function A. In this case 
we can calculate the asymptotic behavior of the Toeplitz 
determinants using Widom's theorem . This theorem 
states that for N± — ► 00 and for m r > — 1: 



detD -> exp{c iVi} N( 



(16) 



Widom's theorem cannot be applied to A, since 
\- 1 / 2 (9) ~ H r (2 - 2cos(<9 - a r ))- mr/2 involves nega- 
tive exponents. We thus employ the alternative expres- 
sion HI containing the matrices D and F. Since the 



elements of D and F can only be obtained by the Fouri- 
ertransform Q if their dimension is at most (N + l)/2, 
there is only one particular decomposition which we can 
consider, namely N± = (N — l)/2 and N 2 = (N + l)/2. 
For the same reason it is not possible to apply Widom's 
theorem to V 1 / 2 as a whole, det V 1 ! 2 can however eas- 
ily be calculated directly from the discrete eigenvalues 
©■ After a lengthy but straight forward calculation we 
eventually obtain the following expression for the mutual 
information with Ni = (N - l)/2 and N 2 = (N + l)/2 



In N + const. 



(17) 



Thus a singular spectral function X 1 ^ 2 {9) in the case of 
half/half partitioning leads to a lower bound to the en- 
tropy that grows logarithmically with the number of os- 
cillators stating a break-down of the area law of entangle- 
ment. As shown above a singular spectral function also 
implies a diverging spatial correlation length, defining a 
critical system. 

The above discussion can be extended to d dimen- 
sions. In this case one would consider the entropy 
S of a hypercube of oscillators with dimensions N\ x 
N2 x • • • x Nd- Since we are interested in the ther- 
modynamic limit we can again assume Ni < (N + 
l)/2. In this case the matrices A and D are Toeplitz 
matrices with respect to each spatial direction and 
their elements A kll k 2 ,...,kd can be obtained from the 
square root of the <i-dimensional function A(#i, . . . , 9d) = 

is regular, the (i-dimensional Szego theorem holds |13| . 
which asserts that the Toeplitz determinant of dimension 
n\ x ni x ■ • • x rid has the asymptotic form 



det D — > exp< cqU\ 



a 

|com---nd + E 



"1 



C 



(18) 



where c = ^ ■ ■ ■ J^d6 d \n( y X 1 / 2 (9 1 , . . . , 9 d )) , 

and the Cj are some constants, whose explicit form 
is of no interest here. We see that under the above 
conditions for the (i-dimensional characteristic function 
A 1 / 2 (^i, . . . , 9d) the entropy has the lower bound 



*>E 



nin 2 ■■■rid 



C, 



,d-l 



(19) 



which is again proportional to the surface area. We note 
that the lower bound (|19|) to the entropy given by the 
multi-dimensional Szego theorem is more general than 
the estimates given in Q and 0, which are restricted 
to nearest neigbor interactions. From the exponential 
bound to the matrix elements of V" 1 ^ 2 one can also find 
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an upper bound to the entropy using eq.JSJ) 



-1/2 



n x n 2 ■ ■ ■ n d 



11 i 



(20) 



Eas. l|19[l and 120() establish an area law for arbitrary di- 
mensions in the case of a regular spectral function. 

In order to obtain a lower bound to the entropy for 
a singular spectral function in more than one dimension 
and to show a corresponding break-down of the entangle- 
ment area law, one would need a multi-dimensional gen- 
eralization of Widoms theorem [l4|. Although no such 
generalization is known to us, there is strong evidence 
for a break down of the area law in higher dimensions. 
First of all for an interaction matrix that is separabel 
in the d dimensions, i.e. whose elements can be written 
as products Vi 1 j 1 Vi 2 j 2 . . . Vi d j d , the ID discussion can 
straight-forwardly extended to d dimensions. Secondly 
Widom has given a generalization of his matrix theorem 
to operator functions f(A) on R d 0]. The proof given 
in [15| makes however use of strong conditions on / that 
are not fulfilled for the case we are interested here. 

To illustrate validity and break-down of the area theo- 
rem let us consider the Hamiltonian H = + 
5 Y^i=i ( — 2r7Qi + q i+ i + qi^i f with periodic boundary 
conditions. The square root of the spectral func- 
tion reads in this case X l / 2 (0) = 1277 — 20086*1. For 
7] > 1, A 1 / 2 is regular and the correlation length is fi- 
nite. For 77 < 1, A 1 / 2 (0) can be written as X 1 / 2 (9) = 
(2-2cos(6» + o )) 1/2 (2-2cos(6»-6' o )) 1/2 , with 77 = 
cos^o, and thus is singular. In this case the correla- 
tion length is infinite. We have numerically calculated 
the entropy for this system for different values of r\. The 
results are shown in fig^ O ne recognizes an unlimited 
logarithmic growth of S for 77 < 1 and a saturation for 
77 > 1. 

In the present paper we discussed the relation be- 
tween entanglement and criticality in translational in- 
variant harmonic lattice systems with finite-range cou- 
plings. We have shown that upper and lower bounds to 
the entropy of entanglement as well as the correlation 
length are solely determined by the analytic properties 
of the spectral function. If the spectral function is regu- 
lar, the entanglement obeys an area law and the system 
is non-critical. If the spectral function has a singular 
part, the area law breaks down and the system is criti- 
cal. Thus for harmonic lattice systems with translational 
invariant, non-random and finite-range couplings there is 
a one-to-one correspondence between entanglement and 
criticality. We note that some of our results apply also to 
more general couplings. For the estimates of the entropy 
it is sufficient that the number of roots of h(z) on the 
unit circle is finite. This is always fulfilled for banded 
coupling matrices V but also holds under more general 
conditions. For couplings of infinite range, the regular 
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FIG. 1: Entropy as function of partition size for non-critical 
(77 = 1.2,1.6) and critical harmonic chain (77 = 0.2,0.6) for 
the example of the text obtained from numerical calculation 
of eq.JTJ. 



part of the spectral function Ao is no longer a polyno- 

— 1/2 

mial. Thus A may not be smooth anymore and could 
have a singularity in a derivative of some order. In such 
a case the spectral function could be regular, allowing 
for an entanglement area theorem, and at the same time 
the correlation length would be infinite, i.e. the system 
would be critical as in the example of Ref . |IJ . 
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